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Primordial magnetic fields could provide an explanation for the galactic magnetic fields observed 
today, in which case they may also leave interesting signals in the CMB and the small-scale matter 
power spectrum. We discuss how to approximately calculate the important non-linear magnetic 
effects within the guise of linear perturbation theory, and calculate the matter and CMB power 
spectra including the SZ contribution. We then use various cosmological datasets to constrain the 
form of the magnetic field power spectrum. Using solely large-scale CMB data (WMAP7, QUaD 
and ACBAR) we find a 95% CL on the variance of the magnetic field at 1 Mpc of B\ < 6.4 nG. 
When we include SPT data to constrain the SZ effect, we find a revised limit of B\ < 4.1 nG. The 
addition of SDSS Lyman-a data lowers this limit even further, roughly constraining the magnetic 
field to B x < 1.3 nG. 



I. INTRODUCTION 

Magnetic fields are observed at the 10 _6 G level in 
galaxies and clusters, not just in the local Universe, but 
up to a redshift of z ~ 0.7-2 [1 j. Production of mag- 
netic fields within a formed galaxy is extremely difficult, 
and the consensus is that they are amplified from pre- 
galactic seed fields [2 . One interesting possibility is that 
these seed fields were primordial in origin, formed in the 
very early Universe. 

The observations of magnetic fields up to a redshift of 
z ~ 2 are particularly interesting as the fields have ampli- 
tudes (B ~ 10/iG) comparable to those observed locally 
[3] . This evidence seems to disfavour a large dynamo am- 
plification, requiring a larger seed field to produce today's 
magnetic fields. In the extreme case of there being no 
dynamo amplification, adiabatic contraction alone could 
amplify a pre-galactic field of around 1 nG (comoving) to 
the required level. 

A stochastic primordial magnetic field would source 
scalar, vector and tensor modes, giving rise to both tem- 
perature and polarisation perturbations in the Cosmic 
Microwave Background [IH8], and perturbing the matter 
distribution in the local universe [9] . Our aim in this pa- 
per is to see what limits current data places on the level 
of magnetic fields in the early universe, and whether this 
is compatible with a scenario where galactic magnetic 
fields are seeded by primordial magnetic fields. Similar 
work has been carried out in the past [10-12 , finding 
field strengths of several nG to be consistent with cur- 
rent CMB data. 

This work will depend heavily on the results of a pre- 
vious paper [8], and will use the same conventions and 
notation. Where it is necessary to use perturbation the- 
ory we use a gauge invariant notation [T3l [14] similar to 
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the Conformal Newtonian Gauge. In this work we limit 
ourselves to a flat ACDM universe. 



II. NON-LINEAR FIELD EVOLUTION 

We will consider a stochastic magnetic field B l (x^ , r) 
generated by some mechanism in the very early Universe. 
As for all the periods of interest the Universe contains a 
highly ionised plasma, we can use the MHD equations to 
describe the behaviour of the magnetic field. At linear 
order in the magnetic field 



[a 2 ^(x,r)j = 0, 



(1) 



and this motivates the definition of a comoving magnetic 
field, B z (-x,r) = a 2 £?*(x, r) which is time independent 
at linear order. Using this we can write the non-zero 
components of the magnetic energy-momentum tensor as 
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As there is no magnetic field on the background, the per- 
turbations of the stochastic background are manifestly 
gauge invariant. We construct density and anisotropic 
stress perturbations, Ab and n#, defined by 

(3a) 
(3b) 
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where we include the factors of p 7 and p 7 to take ac- 
count of the a -4 factors. As usual the anisotropic stress 
I^B l j can be decomposed into scalar, vector and tensor 
contributions. 

At higher order the comoving magnetic field obeys 

d T B l = [Vx (vxB)+r/V 2 B] i 

= e^ k e klm d 5 {y l B m ) + r l d j d j B i (4) 
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where v is the baryon velocity and 77 is the magnetic 
diffusivity. In a highly conductive medium 77 is negligi- 
ble and so we will set it to zero from here onwards. As 
Q is non-linear the magnetic field evolution cannot be 
treated accurately in the standard linearised Einstein- 
Boltzmann approach. The standard approach is to sepa- 
rate the higher order evolution of the magnetic field into 
two effects that are physically well- motivated, a damping 
due to the radiation viscosity, and the effect of the mag- 
netic Jeans length. For extended discussions, see [T5UT7] . 
and especially [18] to which we follow a similar line of 
reasoning below. 

First let us construct A 2 - 7 , a quadratic combination of 
the magnetic field 
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leading to greater pressure opposing to gravitational col- 
lapse, which leads to a magnetic Jeans length. We can 
see that the relevant terms in (191) and (11) are the first 



terms on the right-hand side of each, to this end we will 
largely neglect the other terms. This leaves 



A B = --A B 9 b , 
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We can use the baryon density evolution equation 
A b = -9 b + 3$ 



(12) 
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to replace #5, and neglect 4> on the small scales of interest, 
finding solutions for the magnetic perturbations 



with the normalisation chosen such that it is conveniently 
close to the magnetic perturbations 
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As = 
Ii% = A k k 5 ij - 3A ij 
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We are interested in the time evolution of A ZJ , and thus 
want to calculate the derivative = ~ — - — -^B^ l d T B^\ 
where the parentheses indicate symmetrisation with re- 
spect to the enclosed indices. It will be useful to re- 
express Q by expanding the Levi-Civita symbols 

d T & = B k (d k v i ) - v k (d k B i ) - frdkv 1 * . (7) 

Using ([7]) we can calculate the time derivative of A ZJ 

A ij = 2(d k v^)A j ^ k - 2(d k v k )A ij - v k (d k A ij ) . (8) 

Using the definitions of ([6| we can split the above into 
equations for each of the perturbations A B and U B . 
Firstly 

A B = -\A B h - v\d k A B ) - i<7 fcz llg (9) 

where we have only included terms up to second order, 
and we have decomposed 
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at linear order. Physically #5 is the divergence of the 
baryon 3-velocity, and is its shear; we will neglect 
the antisymmetric vorticity tensor ujij. The magnetic 
anisotropic stress evolves as 
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The first effect to consider is that as perturbations col- 
lapse the magnetic field lines will also be compressed, 



A B = A^oexp ( -A b 



n2 = n£ exp -a 6 



(14) 



where A^o and 11^,0 are the magnetic perturbations at 
some initial time where the baryon perturbation A5 was 
much smaller than its present value. The linearised equa- 
tion for the evolution of the baryon velocities in Fourier 
space is 
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where r c is defined as r c = l/(an e cr T ), c Syb the baryon 
sound speed, R = 4p 7 /3p&, and the last term is the 
magnetic interaction, the Lorentz force. To include the 
magnetic evolution in this, we need to somehow linearise 
(14), which we do by expanding up to second order, and 



replace the magnetic perturbation with its expectation 
(smoothed at the relevant scale) in the higher order terms 
only. In Fourier space this leaves 



A B = A B , + ^(A B ,o) fc A 6 , 



(16) 
(17) 



The equation for n# has only the lowest order term as 
(IIb,o) k = 0. There is some justification for this — as 
A B is always positive, whereas A5 has zero mean, we 
expect there to be a significant coherent term from the 
large-scale magnetic energy density times the baryon per- 
turbation at the scale of interest. 



Inserting (16) into (15) gives 
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where we have defined the Alfven velocity va as 



equations for the two perturbations for r < r* 1 



-2k 2 v\ t c A B , 
-2k 2 v 2 a t c U b . 



(21) 



where p is the density of the conducting fluid. The last 
equality comes from the fact that during matter domina- 
tion p = p b . The overall effect of including the magnetic 
Jeans effect is to appear as the standard evolution equa- 
tion ( 15 ) with a new effective sound speed c 2 —> (? s + \v\. 
This agrees with other approximate treatments in the lit- 
erature up to the factor of | [15j [17] . As this effect is 
non-linear, even when evolving the standard adiabatic 
mode we must add continue to add in the Alfven veloc- 
ity term. 

Radiation free streaming is particularly important for 
magnetic fields, and is the most important source of 
damping on large scales [16 . Prior to free-streaming, 
there are many photon-baryon scatterings per wave- 
length, and the radiation and baryons appear like a 
single tightly-coupled conducting fluid. However, when 
the photons start to free stream they decouple from the 
baryon fluid and this becomes the sole conducting fluid. 
Whilst the baryons are no longer tightly coupled to the 
photons, there is still enough scattering to exert a signif- 
icant drag force on the fluid, and this causes damping on 
propagating magnetic waves. 

When the photons decouple from a particular scale, 
the fluid no longer feels the radiation pressure, and the 
magnetic pressure dominates the baryon pressure for B > 
0.1 nG. In this regime, in addition to Alfven modes, 
magneto- acoustic modes are also significant. 

Properly accounting for the damping in the radiation 
free-streaming regime requires the use of the full Boltz- 
mann system for the photons, a procedure that is compli- 
cated by the necessity of including non-linear magnetic 
effects. Instead we use the prescription of [15], who anal- 
yse the evolution of Alfven modes in the presence of a 
homogeneous radiation drag force. They find that the 
magnetic field on small scales damped approximately as 



B\k) = Bj(k) exp (-k 2 pv\r c dr^j 



(20) 



where va and r c are defined above, and r* is the time of 
last scattering. For magneto-acoustic modes, the damp- 
ing is similar to that of ( 20 ) with the Alfven velocity va 



replaced with the baryon sound speed c s? 5, causing them 
to damp on smaller scales [15j [16] . For simplicity, we will 
treat the field as damping solely on the largest scale, that 
of the Alfven modes. 

To account for this damping in our work, we allow the 
linear perturbations Ab and Ub to evolve. This is in 
contrast to other work in the literature [5j |6j [10] where 
the perturbations are constant for all time and the damp- 
ing is imposed as a cut off on the initial power spectrum. 
This is implemented by constructing effective differential 



For r > r*, we set Ab and Hb to zero, as after recombi- 
nation there is less than one scattering per Hubble time 
and the damping is negligible. To ensure that energy 
conservation is maintained we presume that this is trans- 
ferred equally into the photon-baryon fluid. We modify 
the baryon density equation to 
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and the photon density equation to 
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where we have added the final term in both of the 
above equations. Practically, ensuring energy conserva- 
tion makes little difference to the results. 



III. MAGNETIC MATTER POWER SPECTRUM 

Using the above modifications to the evolution equa- 
tions we are able to calculate matter power spectra in- 
cluding approximate treatments of the important non- 
linear effects. The remaining input is the statistics of 
the initial magnetic perturbations Ab and n#. We as- 
sume that at some early time the comoving magnetic field 
power spectrum is described by a power spectrum 



B i (k)B*(k')) = (2nf5(k 
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P B {k), (24) 



where P ? , 



kikj, and we will not consider helicity. 



The magnetic field power spectrum is 
P B (k) = Ak nB . 



(25) 



We will use the spectral index n# as defined, but it is 
conventional to give the amplitude of the spectrum in 
terms of the variance B 2 of the magnetic field strength 
at a scale A (we choose A = 1 Mpc to agree with the bulk 
of the literature). This gives 



A 



{2n) nB+b B 2 x 



2T (^±3) kl B 



+3 



(26) 



1 Though we base our treatment on that of 15 , it is worth not- 
ing t hat appro xim ately the same result for Ab can be der ived 
from ( |12| ) and (|15|, using the same approximation as in ( |16| ) and 
those within |15) ( not ably, homogeneous baryons ~ and 
drag 6 1 ~ 0, and terminal velocity A# ~ 0) in the overdamped 
regime. 
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FIG. 1: The additional non- linear effects modify the magnetic 
power spectrum. We plot: the linear only behaviour; the 
diffusion damping of the magnetic field; the magnetic pressure 
support; and both effects combined. For comparison we also 
include the primary adiabatic mode. 

To calculate the two required power spectra P& B (k), 
Pu B (k) and the cross spectrum Pa b u b (^) requires con- 
volutions of the magnetic power spectrum. We use the 
results of a previous paper [8] , which allows us to reduce 
the convolution into a dimensionless integral depending 
only on n^. We numerically evaluate these at a large 
number of values such that we can accurately interpolate 
to find power spectra at an arbitrary n#. It is important 
to note that as we are including the damping and Jeans' 
effects in the evolution of the magnetic perturbations we 
can safely avoid imposing a cutoff in the initial power 
spectrum as used in other work [5j [12] . 

In Fig.[T]we show the consequence of the two non-linear 
effects on the matter power spectrum. The first thing to 
note is that the linear theory magnetic power spectrum 
grows as P(k) cx k. On sub-horizon scales during matter 
domination, we can combine the evolution equations for 
baryons and dark matter to give an equation for the total 
matter perturbation A m = RbA^ + R C A C 

A m + UA m - \u 2 A m = ^k 2 L B , (27) 

where the Lorentz force Lb = | (|n# — A#), and we 
have also neglected pressure terms in c 2 s b . In these limits 
the equation has a simple solution for the magnetic mode 

A m = R h L B P -^p\ {knf [1 (^) 2 + • • • 1 (28) 

p h [Ti) LIOVt;/ J 

where we have included only the leading order term, and 
Ti is the time when the mode starts to grow significantly. 
This is the time that the baryon perturbation decouples 
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FIG. 2: The full matter power spectrum with all non-linear ef- 
fects at a variety of magnetic field strengths between B\ — 1- 
5nG. We have held the magnetic spectral index constant at 
riB = —2.9. The power spectrum amplitude increases strongly 
with that of the magnetic field, however the two scales asso- 
ciated with the diffusion damping and magnetic Jeans length 
also increase rapidly. 



from the photons. Provided this time is similar across a 
range of scales (correct for larger scales), we can expect 
Pm{k) cx k 4 PL B (k); for nearly scale invariant magnetic 
field spectra this gives P m (k) cx k. 

The effect of including the magnetic pressure increases 
the Jeans length, and thus causes smaller scales to os- 
cillate. In matter domination, the comoving magnetic 
Jeans wavenumber is constant and thus larger k grow 
very little after recombination. On even larger scales the 
diffusion damping of the magnetic fields causes the source 
for the later growth of the magnetic fields to be exponen- 
tially suppressed. This leads to much slower growth in 
the matter perturbations, and an effective cutoff beyond 
which there are no significant perturbations sourced. The 
scales at which these effects start roughly agree with the 
estimates of [I5j [TTl [19] , where both are expected to scale 
like k c oc B~ (n+5)/2 . 

In Fig. [2] and Fig. |3j we plot the effect on the matter 
power spectrum of the amplitude and tilt of the mag- 
netic spectrum. Fig. [2] shows that the relative contrast 
between the peak magnetic and primary contributions 
stay roughly constant as the amplitude B\ is varied. 
This occurs as a consequence of the fact that the matter 
power spectrum at the damping cutoff is cx B\k ci whilst 
the primary spectrum cx k~ s . However as the damping 
wavenumber k c cx B^ 1 for nearly scale invariant spectra, 
the ratio between the two is constant. 
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FIG. 3: The full matter power spectrum with all no n- linear 
effects at a variety of magnetic spectral indices between nu — 
—2.9 to —1.8. The amplitude is fixed at B\ = 1 nG. 



IV. RESULTS 

We have used versions of CAMB [20] and Cos- 
moMC [21] modified to generate the magnetic con- 
tributions to both the CMB angular power spectrum 
and the matter power spectrum 2 . We sample with a 
flat prior in the standard set of cosmological parame- 
ters Q c h 2 , r, n s , log (10 10 A S )}, assuming a flat 
ACDM universe. To describe the magnetic fields we sup- 
plement this with three parameters, the magnetic power 
spectrum amplitude B\ and spectral index n#. We also 
need to specify the production time of the magnetic fields 
tb to calculate the passive modes. For this we use a proxy 
tb = l°gio( r ^ / t b) ? where r v is the time of neutrino de- 
coupling. 

For priors on the magnetic quantities we choose a flat 
prior in B\ < 10 nG. For we use a flat prior with 
bounds 6 < tb < 17 corresponding to magnetic field pro- 
duction between the electroweak phase transition, and re- 
heating at the GUT scale. Finally for the spectral index 
we use a flat prior —2.95 < n# < —1.6, the lower bound 
of this comes from the fact that there is an infrared di- 
vergence for riB < —3. The upper bound is from the 
fact that we are primarily interested in primordial fields 
produced in the early Universe (prior to nucleosynthe- 
sis), the work of [22] suggests that larger spectral indice s 
are ruled out by nucleosynthesis constraints (see § IV D). 



Other than the upper bound on the spectral index, our 
results should be valid for fields produced at later epochs. 



A. CMB only 

Primordial magnetic fields contribute both passive and 
compensated scalar and tensor modes, as well as a com- 
pensated vector mode. However, the significant contri- 
butions are from the passive tensor mode, and the com- 
pensated vector mode, which are important on large and 
small scales respectively [8 . For comparison to CMB 
data (for scales I < 3000) we need only need calculate 
these two. 

We use the recent WMAP 7 year release [23 along with 
the final data from QUaD [24 and ACBAR [25] for I < 
3000. We also use the Hubble Key Project data [26], BBN 
data, Union Supernova dataset [27] and BAO data from 
SDSS DR7 [28]. To account for the Sunyaev-Zel'dovich 
effect at high multipoles we adopt the standard treatment 
and use the WMAP template [23] . This requires an extra 
parameter for its normalisation Asz which we treat as a 
nuisance parameter and marginalise over. 

The resulting constraints are shown in Fig. [4j which 
are in broad agreement with those of Ref. [12 . We have 
included only the magnetic amplitude B\ and spectral 
index n#, marginalising over the other parameters. The 
standard cosmological parameters are not shown as they 
are in agreement with their values in a universe with no 
primordial magnetic fields. We have also chosen not to 
include r# as it is unconstrained by the data (see Fig.[5|. 
This is a manifestation of the fact that the tensor pas- 
sive mode can be only significant at large magnetic am- 
plitudes, and these are already excluded by the magnetic 
vector mode. From this we calculate a 95% Confidence 
Limit of B x < 6.4 nG. 

The constraints on the magnetic parameter space are 
unchanged when adding in large-scale matter power data 
sets such as 2dF [29] and SDSS LRG data from DR4 [30] 3 , 
though the cosmological parameters change as expected. 
The reason can be seen from Fig. [2] as the magnetic con- 
tributions are only significant for k > 1/iMpc -1 , and the 
galaxy redshift surveys probe only up to k ~ 0.2/iMpc -1 . 
We will include this data for the remainder of this paper. 



B. Sunyaev — Zel'dovich Effect 

Recent work [19] has suggested that the Sunyaev- 
Zel'dovich (SZ) effect may be able to give tight con- 
straints on the magnetic power spectrum. As the mag- 
netic fields cause an increase in the small scale matter 
power spectrum, this gives rise to a large number of small 
mass halos and thus a rise in the SZ angular power spec- 
trum. The recent release of data from the South Pole 



3 We do not use the latest SDSS DR7 data [3j due to complica- 
2 The modified version of CAMB and the adaptation of CosmoMC tions modifying its likelihood calculation to include the magnetic 



to use it will be available from http://camb.info/jrs field matter power spectrum. 
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FIG. 4: The CMB only (WMAP7 + QUAD + ACBAR) con- 
straints on the magnetic field amplitude and spectral index. 
The top left plot is the marginalised constraint solely on the 
magnetic field, and the bottom right is the marginal con- 
straint on its spectral index. The bottom left plot is the joint 
constraint (marginalised over the other cosmological parame- 
ters), with contours for one and two sigma errors. 

Telescope (SPT) [32j, provides CMB observations up to 
/ ~ 10 4 that we will compare to. 

We follow the general method of [19] , using the pre- 
scription of [33 to calculate the thermal SZ angular 
power spectrum from the matter power spectrum, which 
we have calculated using our modified version of CAMB. 
We give the details of this calculation in Appendix [A[ 
In addition to the SZ contribution the magnetic vector 
mode gives a significant addition to the small scale CMB 
power, which must be included. We do not include the ki- 
netic SZ spectrum which is sub-dominant to the thermal 
spectrum. We also note that the non-gaussian statistics 
of the SZ effect can significantly increase its intrinsic vari- 
ance [34.; we do not account for this within this paper. 

As well as the datasets for the CMB-only constraints, 
and the SPT data, we have also included the large scale 
matter power data. When generating the SZ power spec- 
trum we must include all contributions, both the adia- 
batic and magnetic contributions, as such we do not add 
in an SZ template, and thus Asz is no longer sampled 
over. The SPT data contains a further small scale contri- 
bution from the unresolved point sources, both from star- 
forming and radio galaxies (see Fig. [6]). The exact value 
must be fit from the bandpowers (in the analysis of [32] it 
is poorly constrained with Cf s = (6.2 ±6.4) x l(T 7 /iK 2 ). 
We expect that neglecting this contribution will have re- 
sulted in a small increase to our upper limits. 

The magnetic contribution to the SZ power spectrum 
from our calculations is smaller than that of [19 j. Dif- 
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FIG. 5: The CMB-only marginal constraint on the produc- 
tion time ratio tb — log 10 (t v /tb) (bottom left), along with 
the correlations with the magnetic field amplitude (top left) 
and spectral index (bottom right). Even in the CMB-only 
constraint tb is largely unconstrained and has minimal cor- 
relation with the other parameters, being only important in 
the case of very large magnetic fields, and a very red spectral 
index. 

ferences between our calculations, such as in finding the 
amplitude of the matter perturbations and treating the 
non-linear effects, make it difficult to give a single reason 
for this. 

Figure [7] shows the marginalised probability distribu- 
tions for the magnetic parameters (again we do not plot 
tb)- The addition of the SPT data has excluded much 
of the parameter space that was allowed when includ- 
ing primary CMB effects only, especially the region with 
large magnetic field and blue spectral index, which gives 
the most dramatic change in the matter power spectrum 
(see Fig. [3|. This gives a large decrease in the 95% CL 
of the magnetic amplitude to B\ < 4.1 nG. 



C. Lyman Alpha Data 

It is clear from both Fig. [2] and Fig. [3] that the small 
scale matter power spectrum is significantly affected by 
the presence of a primordial magnetic field. Unfortu- 
nately galaxy redshift surveys such as SDSS probe only 
as far as k < 0.2 h Mpc -1 , too large scales to be af- 
fected by the magnetic fields. However, observations of 
the Lyman-Qf flux power spectrum probe the matter den- 
sity power spectrum to scales as small as k ~ 5 /iMpc -1 , 
and may be able to give a much more powerful constraint. 
Unfortunately there is no simple analytic mapping from 
the observations to the matter spectrum, so cosmological 
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FIG. 6: The angular power spectrum of the CMB and the 
SZ contributions to it. We plot the SZ contribution from the 
primary adiabatic mode only, and in combination with four 
magnetic field strengths B\ = 2-5 nG, with a magnetic spec- 
tral index of jib = —2.4. We also plot the bandpowers from 
SPT, and the estimate of the residual point source contribu- 
tion (labelled PS). Both the SZ contribution and the SPT 
data are at 150 GHz. 




FIG. 8: Constraints on the magnetic field amplitude and 
spectral index using CMB and matter power data, includ- 
ing Lyman-a data. The Lyman-a data, which probes the 
small scale matter distribution, has dramatically reduced the 
allowable parameter space, to a range of small amplitudes 
(B\ < 1.5 nG) roughly independent of spectral index ub- 
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FIG. 7: The CMB and SZ constraints on the magnetic field 
amplitude and spectral index. The individual plots are equiv- 
alent to those in Fig. [4] Calculating the magnetic SZ contri- 
butions with the addition of SPT data has restricted the pa- 
rameter space to lower amplitudes and redder spectral indices 
than the large scale CMB data. 



constraints need to be obtained by comparison to simula- 
tions [35] . A fully consistent analysis including magnetic 
fields is beyond the scope of this paper; instead we use a 
very rough simple approximation to the likelihood. 

We use the standard Lyman-a module in CosmoMC 
which is based on the work of [35]. This finds an effec- 
tive amplitude and spectral index about a pivot scale of 
k = 0.009 skm -1 (roughly k ~ 1 h Mpc -1 ). It calcu- 
lates a likelihood from these by interpolating between a 
set of simulations compared to SDSS quasar data. Due 
to the large difference between our power spectra (when 
including magnetic effects) and those of ACDM we can- 
not expect the likelihoods to be very accurate, especially 
at large magnetic amplitudes and blue spectral indices, 
where the effect is greatest. However for the pivot scale, 
and range of scales probed by the SDSS spectra used 
k < 0.02 skm -1 , the magnetic contribution is generally 
small compared to that of the primary adiabatic mode. 
In this light our results should be viewed as an approxi- 
mation to the constraints that a more sophisticated like- 
lihood approach would achieve. 

The magnetic parameter space is significantly con- 
stricted by the use of the Lyman-a data (see Fig. [8|, 
with the allowed region for B\ being largely indepen- 
dent of the spectral index ub- As the Lyman-a pivot 
scale (k ~ 1 h Mpc -1 ) coincides with the scale at which 
the magnetic matter power spectrum amplitudes are very 
similar across a large range of n# this is to be expected. 
Overall this results in a 95% CL of B x < 1.3 nG. 
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FIG. 9: Constraints on the magnetic parameter space from 
Lyman-Q! data including the nucleosynthesis constraints of 
[22] . The solid and dashed orange lines give the upper bounds 
in amplitude of a magnetic field produced at GUT scale infla- 
tion and the electroweak phase transition respectively This 
corresponds to our priors. 



D. Nucleosynthesis Constraints 

Nucleosynthesis places strong constraints on the 
amount of energy density in gravitational waves allowed 
in the Universe, giving a limit of ^gw ^ 1-1 x 10 -6 [36] . 
As the anisotropic stress of an inhomogeneous magnetic 
field sources gravitational waves, this allows the small 
scale magnetic fields to transfer their energy into grav- 
itational waves before being dissipated [22]. The limit 
on (]qw can therefore be used to constrain the mag- 
netic field parameter space. We take Eq. (33) from [22] 
adapted to our conventions, yielding 



B x I nG < 700 h 



[ r(f) 



1/2 



10" 



-(4+r B ) 



(29) 



It is essential to note that this assumes that the magnetic 
field is well described by a single power law across a vast 
range of scales — from the pivot scale at k = IMpc -1 up 
to at least k — 10 10 Mpc _1 , the horizon scale at the elec- 
troweak phase transition. Depending on how the mag- 
netic field is generated this assumption may break down. 

In Fig. [9] we show the probability distribution of the 
magnetic parameters when combining the nucleosynthe- 
sis constraint with the Lyman-a data of § IV C As we 



would expect this reduces the allowed parameter spaces 
to redder spectral indices n# < —2. Because of this the 
limits on the amplitude are slightly enlarged with the 
95% CL becoming B x < 1.6 nG. 



V. CONCLUSION 

As we have seen in this paper, when considering the 
magnetic contributions to the matter power spectrum it 
is essential to treat important small-scale, non-linear ef- 
fects. We have demonstrated a technique for approximat- 
ing the main non-linear effects within linear perturbation 
theory, and have incorporated this into a modified ver- 
sion of CAMB. This gives an alternative to the common 
approach of incorporating these effects directly into the 
initial power spectrum, such as in Ref. [T2] . 

We have used our theoretical predictions to place con- 
straints on the allowable magnetic field amplitude given 
various data sets. Where all limits are at 95% confi- 
dence, using CMB data only we find B\ < 6.4 nG, with 
a redder spectral index favoured. As the presence of a 
stochastic magnetic field gives significant modifications 
to the small scale matter distribution, we also look at the 
constraints when adding two probes of it: the Sunyaev- 
Zel'dovich effect measured by SPT [32], which gives a 
limit of B\ < 4.1 nG; and Lyman-a forest data from 
SDSS [35], which gives a rough constraint of B\ < 1.3nG. 

Whilst the addition of the small-scale matter data gives 
a large reduction in the allowed amplitude of a primor- 
dial magnetic field (certainly when using Lyman-a data), 
this is still roughly consistent with a scenario where cur- 
rent galactic fields are formed solely by adiabatic con- 
traction of primordial fields. Increasing the accuracy of 
measurements at the scales we can currently probe with 
data will provide limited gains: since the magnetic power 
spectrum increases with B^ significantly decreasing the 
errors on matter power spectrum measurements will pro- 
duce less impressive gains in the upper limit on B\ . How- 
ever, if better observations and modelling allow accurate 
comparison to smaller scales there should be an almost 
linear decrease in the upper limit with the smallest scale 
probed. This would seem to provide the best opportunity 
for testing the primordial field hypothesis. 

Our modifications to CAMB to for calculating the 
magnetic matter CMB and matter power spectra, as well 
as the general halo mass function and SZ power spectrum 
will shortly be available from |http: //camb. inf o/jrs. 
The modifications to integrate them into CosmoMC will 
be placed at the same location. 
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Appendix A: Sunyaev — Zel'dovich Effect 

In order to compare to the recent data from the South 
Pole Telescope [32], we need to be able to calculate SZ 
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angular power spectra from linear matter power spectra. 
We use the halo method of Komatsu and Seljak [33 , and 
largely follow details of the calculation in [32 j. We give 
an outline of this below. 

The angular power spectrum is given by 



The prone y 3 D is given by 

Vsd(x) = 1.14 x 10" 4 Mpc _1 
Pb(0) 



k B T b (0) 
8keV 



y(x). (A7) 



clM^§^\y l{ M, Z )f . (Al) 



In the above V c is the comoving volume out to redshift 
z, and g v is the spectral function given by 



LlO 14 M Mp C - 3 . 
The dimensionless function y{x) gives the profile shape 

17/(7-1) 



y(x)= l-5[l-x _1 ln(l + x)] 



9v 



tanh(x/2) 



-4 



(A2) 



where x = hu/ksTcMB- The halo mass function ^ 
is comoving number density of virialized halos at mass 
M. Finally yi is the Fourier transform of the projected 
Compton y profile 



B is a constant given by 
-17-1 



B = 3t]~ 



7 



1 , , ^ 1 

-In 1 + c - 

c v ' 1 + c 



(A8) 



(A9) 



Fitting functions for r] and 7 are derived in [33]. They 
are valid for the range 1 < c < 25: 



Vi(M, z) 



y 3 D(x) sine (lx/l s )x 2 dx , (A3) 7 = 1.137 + 0.0894 In (c/5) - 3.68 x 10" 3 (c - 5) , 



where r s is the scale radius of the profile, l s = d a /r s is its 
angular projection (d a is the angular diameter distance 
to redshift z) and y^D is the three-dimensional Comp- 
ton profile in terms of x = r/r s . The profile ysD is 
determined by the model chosen for the baryon density 
and temperature profile of the halo. To find it we follow 
the details of [33], fixing its form with four assumptions: 
the dark matter density profile is NFW [37]; hydrostatic 
equilibrium between the gas pressure and the halo self- 
gravity; baryon density traces the dark matter density in 
the outer halo; and the gas has a poly tropic equation of 
state oc pj. The results of this are given below, for 
details see [35] . 

Using an NFW dark matter profile, the scale radius 
above r s is the usual NFW definition, r s = r v [ r /c where 
c is the concentration, and r v i r is the virial radius given 

by 



(A10) 

77 = 2.235 + 0.202 (c - 5) - 1.16 x HT 3 (c - 5) 2 . (All) 
The central gas density Pb(0) is 

M 



p b (0) = 7.96 x 10 12 M o Mpc" 



Q b \ y{x) 



-1/7 



-]n(i + c)-— !— 

c K } 1+c 



/ r v ir V 

10 14 M o y V M P C / 
1 



n m j (1 + c) 2 

and the central temperature T^(0) is 

M 



T h (0) = 0.880 keV rj 



( r vir \ 1 
10 14 M o y V M P C / 



(A12) 



(A13) 



M 



47T A c p c 



1/3 



(A4) 



where the virialisation parameter A c can be calculated 
from the spherical collapse of a top hat perturbation. A 
fitting formula for A c is calculated in [38] 



The final ingredient needed to calculate the SZ power 
spectrum is the mass function. In common with [32j [33] 
we use the Jenkins mass function [40] calculated from 
N-body simulations. As with the Press-Schechter pre- 
scription the key quantity is the smoothed variance cr(R) 
defined by 



\R;z) = J 



d\nkW 2 R {k)V{k;z) . 



(A14) 



A c (» = 18tt 2 - 82ft A (z) + 39f2 A 0) 



(A5) 



which is accurate in the range £7 a < 0.9. 

The concentration parameter c which defines the scale 
radius of the profile can be fitted from simulations. We 
use the relation of 39 which takes the form 



c(M,z) 



7.85 / M 



-0.081 



In our work we choose a window function Wn(k) that is 
a spherical top hat in real space. The mass enclosed in 
this comoving scale is simply M = 47rp m5 o^R 3 /3, and this 
defines an obvious mapping between a mass smoothed 
a(M) and <j(R). In terms of <j(M) the mass function of 
0U] is 



(A6) 



(1 + z) - 71 \M* 
where the pivot mass is fixed to be M* = 2 x 10 12 /i _1 M . 



M 2 dn(M- z) 
p m ,o dM 



0.301 



dine 



d\nM 



exp i 



|0.64-ln<j|) 
(A15) 
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